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ENUMERATION OF THE FACETS OF CUT POLYTOPES OVER 
SOME HIGHLY SYMMETRIC GRAPHS 

MICHEL DEZA AND MATHIEU DUTOUR SIKIRIC 


Abstract. We report here a computation giving the complete list of facets for 
the cut poly topes over several very symmetric graphs with 15 — 30 edges, in¬ 
cluding Ks, ^ 3 , 3 , 3 , Xi ? 4 ; 4 , K 5 , 5 , some other Ki m , K\ i m , Prism 7 , APrismQ, 
Mobius ladder M 14 , Dodecahedron, Heawood and Petersen graphs. 

For K%, it shows that the huge lists of facets of the cut polytope CUTPs 
and cut cone CUTs, given in ED is complete. We also confirm the conjecture 
that any facet of CUTPs is adjacent to a triangle facet. 

The lists of facets for K\ 1 m with (£, m) = (4, 4), (3, 5), (3,4) solve problems 
(see, for example, ES]) in quantum information theory. 


1. Introduction 


Polyhedra, generated by cuts and by finite metrics are central objects of Dis¬ 
crete Mathematics; see, say, PH- In particular, they are tightly connected with 
the well-known NP-hard optimization problems such as the max-cut problem and 
the unconstrained quadratic 0,1 programming problem. To find their (mostly un¬ 
known) facets is one approach to this problem. It is a difficult problem since the 
number of facets grows very fast with the size of the problem. However, when full 
computation is unfeasible, having a partial list of facets can help branch-and-bound 
strategies that are then commonly used. Here we consider the case of cut polytope 
over graphs. 

Given a graph G = (V, E) with n = | Vj, for a vertex subset S C V = {1,..., n}, 
the cut semimetric 5s{G) is a vector (actually, a symmetric {0, l}-matrix) defined 
as 


5s{x,y) 


1 if xy € E and IS 1 fl {x : y}\ = 1 
0 otherwise. 


So, 5s can be seen also as the adjacency matrix of a cut (into S and S) subgraph of 
G. Clearly, <J{i,... )n }_s = 5s■ If G is connected, which will be the case in our work, 
there are exactly 2 n ~ 1 distinct cut semimetrics. The cut polytope CUTP(G) and 
the cut cone CUT(G) are defined as the convex hull of all such semimetrics and 
the positive span of all non-zero ones among them, respectively. Their number of 
vertices, respectively extreme rays is 2 n , respectively 2 n ~ 1 — 1 and their dimension 
is |E|, i.e., the number of edges of G. 

The most interesting and complicated case is CUTP(iC„) and CUT(A n ), denoted 
simply CUTP„, CUT n and called the cut polytope and the cut cone. In fact, CUT„ 
is the set of all n-vertex semimetrics, which embed isometrically into some metric 
space h, and rational-valued elements of CUT n correspond exactly to the n -vertex 
semimetrics, which embed isometrically, up to a scale A £ N, into the path metric of 
some iV-cube ■ It shows importance of this cone in Analysis and Combinatorics. 
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The metric cone MET„ is the set of all semimetrics on n points , i.e., those of 
above functions, which satisfy all triangle inequalities , i.e., dij < dik + dkj■ The 
bounding of MET re by Q) perimeter inequalities djj + d,^ + djk < 2 produces 
the metric polytope METP,,. We have the evident inclusions CUT„ C MET n 
and CUTP n C METP„ with CUT„ = MET n and CUTP„ = METP n only for 
3 < n < 4. Another relaxation of the cut-polytope is the hypermetric polytope 
considered in [16]. In Table[l]we give the number of facets, vertices (or extreme rays) 
and the number of orbits for CUTP n , METP„ with n < 8 and for the corresponding 
cones. 

The symmetry group of a graph G = (V,E) induces symmetry of CUTP(G). 
For any U C {1,..., n}, the map 5s ^ Sjjas also defines a symmetry of CUTP(G). 
Together those form the restricted symmetry group AAes(CUTP(G)), the order 
of which is Aut(G)|. The full symmetry group Aut(CUTP(G)) may be 

larger. For n ± 4, Aut(CUTP(AT n )) = ARes(CUTP(AT„)) but Aut(CUTP(A' 4 )) = 
Aut(A' 4 j 4 ) ([IS]). So, | Aut(CUTP(A' n ))| is 2 n ~ 1 n\ for n ^ 4 and 6 x 2 3 4! for n = 4. 
Let A(G) denote 2 1 _ l y l | Aut(CUT(G))|. One can check the following. 

(1) If G is a complete multipartite graph with ti > 1 parts of size oq for \ <i < r 
and 1 < ai < 02 < • ■ • < a r , then \Aut(G)\ = n)=i 

(2) If G is Prism m , APrism m or Mobius ladder M 2 m , then A(G) = \Aut(G)\ 
and, moreover, A(G) = 4ro = 2|P| if m ^ 4, m > 4 and m > 4, respectively. 

(3) Examples of infinite series of singular , i.e., having |Aut(CUTP(G))| > 
|AAes(CUTP(G))|, graphs G are A^,™, ATip^ for any m > 2 and P m: 
C m for any m > 3. All 4-, 5- and 6 -vertex connected graphs are singular, 
except A'i j 3 , ATi^, AT 5 , A '5 — e, AT 5 — ei — e 2 with disjoint edges ei, e 2 and 
19 (among all 112) 6 -vertex graphs. 

In Table [2] we list information on the cut polytopes of several graphs. Full data 
sets are available from [18] . The computational techniques used are explained in 
Section [2] they can be applied to any polytope having a large symmetry group. 
In Section OH the results for CUTPg are detailed. In Section [4] the results for the 
correlation polytopes K n ^ m are presented. Finally, in Section [5] we explain the 
results for graphs without AT 5 minors. 

2. Computational methods 

The second author has developed over the years an effective computer program 
( [19]) for enumerating facets of polytopes which are symmetric. The technique 
used is adjacency decomposition method , originally introduced in El and applied 
to the Transporting Salesman polytope, the Linear Ordering polytope and the cut 
polytope. The algorithm is detailed in Algorithm [l] and surveyed in [ID] . 

The initial facet of the polytope P is obtained via linear programming. The tests 
of equivalence are done via the GAP functionality of permutation group and their 
implementation of partition backtrack. The problematic aspect is computing the 
facets adjacent to a facet. This is itself a dual description problem for the polytope 
defined by the facet F. 

An interesting problem is the check when 1Z is complete. Of course, if all the 
orbits are treated, i.e., if 1Z = T>, then the computation is complete. However, 
sometimes we can conclude before that: 

Theorem 1. Let G(P) be the skeleton graph of a m-dimensional polytope. 
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Table 1. The number of facets and vertices (or extreme rays) in 
the cut and metric polytopes (or cones) for n < 8 . The enumera¬ 
tion of orbits of facets of CUT„ and CUTP„ for n < 7 was done 
in EiElCH] for n = 5, 6 and 7, respectively. The enumeration 
of orbits of extreme rays of MET 7 was done in [2^. The orbits of 
vertices of METP„ were enumerated in m for n = 7 and in [14] 
for n = 8 . For n < 6 such enumeration is easy. 


p 

n = 3 

n = 4 

n = 5 

n = 6 

n = 7 

n = 8 

CUTP n , e 

4(1) 

8 ( 1 ) 

16(1) 

32(1) 

64(1) 

128(1) 

CUTP n , f 

4(1) 

16(1) 

56(2) 

368(3) 

116,764(11) 

217,093,472(147) 

CUT„, e 

3(1) 

7(2) 

15(2) 

31(3) 

63(3) 

127(4) 

CUT n , f 

3(1) 

12 ( 1 ) 

40(2) 

210(4) 

38,780(36) 

49,604,520(2,169) 

MET ra , e 

3(1) 

7(2) 

25(3) 

296(7) 

55,226(46) 

119,269,588(3,918) 

met„, / 

3(1) 

12 ( 1 ) 

30(1) 

60(1) 

105(1) 

168(1) 

METP„, e 

4(1) 

8 ( 1 ) 

32(2) 

554(3) 

275,840(13) 

1,550,825,600(533) 

METP n , / 

4(1) 

16(1) 

40(1) 

80(1) 

140(1) 

224(1) 



|^ 7 | 

Z \ 


a) A cube and the remaining 
connected graph obtained from 
Theorem [T] (i) 


b) A cube and the remaining 
connected graph obtained from 
Theorem [T] (ii) 


Figure 1. Illustration of connectivity results of Theorem [l] 


( 1 ) m) g(p) is at least m-connected. 

(ii) If we remove all the edges contained in a given face F, then the remaining 
graph is still connected. 

Hence, if the total set of facets, equivalent to a facet in 7 Z\T>, has size at most 
m — 1 or contains a common vertex, then we can conclude that 1Z is complete. 
The requirements for applying any of those two criterions are rather severe, but 
they are easy to check and, when applicable, the benefits are large. The geometry 
underlying Theorem [I] is illustrated in Figure [Q 

The adjacency decomposition method works well when the polytope is sym¬ 
metric but it still relies on computing a dual description. This is easy when the 
incidence of the facet is low but become more and more problematic for facets of 
large incidence. However, quite often high incidence facets also have large symme¬ 
try groups. Therefore, a natural extension of the technique is to apply the method 
recursively and so, obtain the recursive adjacency decomposition method , which is 
again surveyed in ra¬ 
in order to work correctly, the method requires several ingredients. One is the 
ability to compute easily automorphism group of the polytope, see [9] for details. 
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Table 2. The number of facets of the cut polytopes 
CUTP(G) of some graphs. In 4-th column, A(G ) denote 
2 1 -l y l| Aut(CUTP(G))|. The case A(G) > | Aut(G)| is indicated 
by* 


G = (V, E) 


\E\ 

A(G) 

Number of facets (orbits) 

K 8 

8 

28 

8! 

217,093,472(147) 

K3,3,3 

9 

27 


624,406,788(2015) 

ifl, 4,4 

9 

24 

2(4!) 2 

36,391,264(175) 

if 1,3,5 

9 

23 

3!5! 

71,340(7) 

K\ 34 

8 

19 

3!4! 

12,480(6) 

AA.3,3 

7 

15 

2(3!) 2 

684(3) 

ifl,1,3,3 

8 

21 

4(3!) 2 

432,552(50) 

ifi,2,m, rn> 2 

m + 3 

3to + 2 

1 Aut(ATi, 2 ,m)| 

8 m + 8 (”)( 2 ) 

P-5,5 

10 

25 

2(5 r 

16,482,678,610(1,282) 

iu, 7 

11 

28 

4!7! 

271,596,584(15) 

K 4,6 

10 

24 

4!6! 

23,179,008(12) 

ifl,5 

9 

20 

4!5! 

983,560(8) 

if4,4 

8 

16 

2(4!) 2 

27, 968(4) 

if 3 ,m,m > 3 

m + 3 

3 m 

| Aut(A' 3 , m )| 

6 m + 24(“)(2) 

A' 2 ,m, m > 3 

m + 2 

2m 

2 m_1 m! Aut(A' 2 ,m) * 

4m 2 (l) 

A'i,m, m > 2 

m + 1 

m 

ml 

2 m(l) 

Km-\- 2 Km-> TTl ^ 2 

m + 2 

2m + 1 

2 m_1 m! Aut(AT lj i jm )| * 

4m(l) 

Km+ 3 A m , 777/ 2 

m + 3 

3m + 3 

3!m! 

4+ 12m(2) 

Km-\- 4 A m , 777/ ^ 2 

m + 4 

4m + 6 

4!m! 

8 ( 8 m 2 — 3m + 2)(4) 

if 8 - ^3 

8 

25 

360 

2,685,152(82) 

/\ 7 - if 2 

7 

20 

240 

31,400(17) 

Dodecahedron 

20 

30 

120 

23,804(5) 

Icosahedron 

12 

30 

120 

1,552(4) 

Cube 

8 

12 

48 

200(3) 

Cuboctahedron 

12 

24 

48 

1,360(5) 

Tr. Tetrahedron 

12 

18 

24 

540(4) 

APrisme 

12 

24 

24 

2,032(5) 

Prism-j 

14 

21 

28 

7,394(6) 

Pyr(Prism> 5 ) 

11 

25 

20 

208,132(22) 

Pyr(APrismi) 

9 

24 

16 

389,104(17) 

Mobius ladder M 14 

14 

21 

28 

369,506(9) 

Heawood graph 

14 

21 

336 

5,361,194(9) 

Petersen graph 

10 

15 

120 

3,614(4) 


Another is good heuristics for deciding when to apply the method recursively or 
not. Yet another is a storing system for keeping dual description that may be 
reused, and this again depends on some heuristics. The method has been applied 
successfully on numerous problems [201 EH [ 22 ] and here. 

The framework, that we have defined above, can be applied, in order to sample 
facets of a polytope, say, P. A workable idea is to use linear programming and then 
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Data: Polytope P and a group G 

Result: Set 1Z of all inequivalent representative of facets of P for G 
F t— facet of P. 

TZ^{F}. 

V <- 0 

while 1Z is not complete do 
F a facet in 1Z\T>. 

DuDU {F}. 

F <— facets of P adjacent to F. 
for F £ J- do 
test <— true 

if F is not equivalent to a facet in 1Z then 
| TZ^1ZU{F}. 

end 

end 

end 

Algorithm 1: The adjacency decomposition method 


get some facets of P. But doing so, we overwhelmingly get facets of high incidence, 
while we may be interested in obtaining facets of low incidence. One can adapt 
the adjacency decomposition method to do such a sampling. Let us call two facets 
equivalent if their incidence is the same. By doing so, we remove the combinatorial 
explosion, which is the main difficulty of such dual-description problem. At the 
end, we get a number of orbits of facets of different incidence, which give an idea 
of the complexity of the polytope. 

3. The facets of CUTP 8 

In El. the adjacency decomposition method was introduced and was applied 
to the Transporting Salesman polytope, the Linear Ordering polytope and the cut 
polytope. For the cut polytope CUTPg, the authors found 147 orbits, consisting 
of 217,093,472 facets, but this list was potentially incomplete, since their method 
did not treat the triangle, pentagonal and 7-gonal inequalities (defined in Section 
3) at that time. Therefore, they only prove that the number of orbits is at least 
147. The enumeration is used in 13® for work in quantum mechanics. 

Sometimes (13 CP), this is incorrectly understood and it is reported that the 
number of orbits is exactly 147 with El as a reference. Here we show that Christof- 
Reinelt’s list is complete. We had to treat the 3 remaining orbits of facets since we 
could not apply Theorem [T] However, we could apply the theorem in deeper levels 
of the recursive adjacency method and this made the enumeration faster. 

Conjecture 1. 7 |12l ll7\ ) Any facet of CUTP ra is adjacent to a triangle inequality 
facet. 

The conjecture was checked for n < 7; here we confirm it for n = 8. 

Looking for a counterexample to this conjecture, we applied our sampling frame¬ 
work to CUTP n for n = 10, 11 and 12. We got initial facets of low incidence and 
then we complemented this with random walks in the set of all facets. This allowed 
us to find many simplicial facets (more than 10,000 for each) of these CUTP n but 
all of them were adjacent to at least one triangle inequality facet. 
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As a corollary to the confirmation of above conjecture for n = 8, we obtain that 
the ridge graphs of CUTPs and CUTg have diameter 4. In fact, the subgraphs, 
formed by the triangle facets, have diameter 2. 


4. Correlation polytopes of K n ^ m 

In quantum physics and quantum information theory, Bell inequalities, involving 
joint probabilities of two probabilistic events, are exactly inequalities valid for the 
correlation polytope CORP(G) (called also Boolean quadric polytope BQP(G )) of 
a graph, say, G. In particular, CORP(A' rlj m) is seen in quantum theory as the set 
of possible results of a series of Bell experiments with a non-entangled (separable) 
quantum state shared by two distant parties, where one party (Alice) has n choices 
of possible two-valued measurements and the other party (Bob) has m choices. 
Here, a valid inequality of CORP(AT rajT7l ) is called a Bell inequality and if facet 
inducing, a tight Bell inequality. This polytope is equivalent (linearly isomorphic 
via the covariance map) to the cut polytope of AA,n,m |171 Section 5.2]. Similarly, 
CUTP(A'i inim) j) represents three-party Bell inequalities. 

The symmetry group of CORP(AT n)ra ) and CUTP(/\ i t n,m) is of order 2 1+n+m n!m!. 

We computed the facets of CORP(AT„ im ) having (n, m) = (2,m) with m < 6, 
(3,4), (3, 5), (4,4) and confirmed known enumerations for (n, m) = (2, 2), (3,3). In 
fact, the page [23] collects the progress on finding Bell inequalities. The case of 
CORP(R njTl ) is called there (2,n, 2 )-setting. The cases n = 2, 3 were settled in [ 23 ] 
and [26], respectively. For n = 4, partial lists of facets were known; our 175 orbits 
of 36,391,264 facets of CORP(AT4 j 4) finalize this case. 

In contrast to the Bell’s inequalities, which probe entanglement between spatially- 
separated systems, the Leggett-Garg inequalities test the correlations of a single 
system measured at different times. The polytope, defined by those inequalities for 
n observables, is, actually (0), the cut polytope CUTP n . 


5. Cut polytopes of graphs without minor K 5 

Given a graph G = (V,E), an edge (iqi^) G E defines an edge inequality 
x(vi,V 2 ) > 0. Similarly, an s-cycle (vi,... ,v s ) of G with s > 3 defines a cycle 
inequality 

s-l 

'^2,x{v i ,v i+1 ) - x[vi,v a ) > 0. 

i=l 

These inequalities and their switching define valid inequalities on CUTP(G). There 
are 2|A| edge inequalities and their incidence is 2^ I -2 . Each s-cycle of G gives 2 s-1 
s-cycle faces, which are of incidence s21 v l~ s . 

In fact, mm and j2?b gives, respectively, that: 

(i) only chordless s-cycles and edges not containing in a 3-cycle produces facets; 

(ii) no other facets exist if and only if G is a A'5-minor-free graph. 

All JC5-minor-free graphs in Table [2] are Aii,2,m; 7\2,m; Km+i — K m with m > 
2,1 < i < 3 and the skeletons of regular and semiregular polyhedra (Dodecahedron, 
Icosahedron, Cube, Cuboctahedron, truncated Tetrahedron, APrismQ^Prism-j). 
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